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Abstract. In this paper, we prove a new identity for values of the Hurwitz zeta 
function which contains as particular cases Koecher's identity for odd zeta values, 
the Bailcy-Borwcin-Bradley identity for even zeta values and many other interesting 
£-H ! formulas related to values of the Hurwitz zeta function. We also get an extension of 

the bivariate identity of Cohen to values of the Hurwitz zeta function. The main tool 
we use here is a construction of new Markov- WZ pairs. As application of our results, 
we prove several conjectures on supercongruences proposed by J. Guillera, W. Zudilin, 



d ■ and Z.-W. S 



un. 



1. Introduction 

t> . The Riemann zeta function for Re(s) > 1 is defined by the series 

in 

m 



where 



F 



cw=E^=-' , .G::::;l| 1 )' 

n=l \ ' ' / 

(li) n ■ ■ ■ (a p ) n 



CL\ , . . . , ftp 

h,...,b q 



oo 



^ (6l)n " • • (b q ) n n\ 

is the generalized hypergeometric function and (a) n is the shifted factorial defined by 
(a) n = ci(ft + 1) • • • (a + n — 1), n > 1, and (a)o = 1. 



In 1978, R. Apery used the faster convergent series for £(3) 



(i) co) = ^E ' 



k- 



to derive the irrationality of this number [21]. The series ([T]) first obtained by A. A. Mar- 
kov [TJj in 1890 converges exponentially faster than the original series for C(3), since 
by Stirling's formula, 

(&-»+oo). 



1991 Mathematics Subject Classification. 33F10, 05A30, 11B65, 05A19, 33D15, 11M06. 
Key words and phrases. Zeta values, Apery-like series, generating function, convergence accelera- 
tion, Markov- Wilf-Zeilberger method, Markov- WZ pair, supercongruences. 

1 This research was in part supported by a grant from IPM (No. 89110024). 

2 This research was in part supported by a grant from IPM (No. 89110025). 

1 



A general formula giving analogous Apery-like series for all ((2n + 3), n > 0, was 
proved by Koecher [T3] (and independently in an expanded form by Leshchiner [TB]). 
For \a\ < 1, it reads 

n=0 fe=l v y fe=l V fc ) m=l v 7 

Expanding the right-hand side of ([2]) by powers of a 2 and comparing coefficients of a 2n 
on both sides leads to the Apery-like series for ((2n + 3) 



(3) C (2„ + 3) = I ± ^ (-D"ef (*) + 2£E jtk_(-xr^(t), 

fc=l V fc / j=l k=l \k ) 

where for positive integers r, s, 

fc-i 

e«(*) := [f][](l+ri) = E Uih---jr)- s , 

3=1 l<jl<32<-<jr<k-l 

and [t r ] means the coefficient of t r . In particular, substituting n = in ([3]) recovers 
Markov's formula (JT|) and setting n — 1, 2 gives the following two formulas: 

(4) c(5) = 2 Sw^S<rS? 

(5) c(7) = 2 Sw _2 Sw^ +i SwS^S?' 

respectively. In 1996, inspired by this result, J. Borwein and D. Bradley [4] applied 
extensive computer searches on the base of integer relations algorithms looking for 
additional zeta identities of this sort. This led to the discovery of the new identity 

which is simpler than Koecher's formula for ((7), and similar identities for ((9), C(H), 
C(13), etc. This allowed them to conjecture that certain of these identities, namely 
those for C(4n + 3) are given by the following generating function formula [3]: 

n=0 fe=l fc=l lfcA A ' " I m=l v 7 

The validity of (J7|) was proved later by G. Almkvist and A. Granville [Tj in 1999. 
Expanding the right-hand side of fl7|) in powers of a 4 gives the following Apery-like 
series for ((in + 3) j3]: 

n oo (_-,\k+l n -3 

( 8 ) ^+ 3 ) ^EE^E 4r O)4 4) w> 

j=0 fc=l V fc / r=0 



fe+1 n oo / 1 s 



where 

fc-i 

hi s \k) [f]H(i - r s t)-\ 

3=1 

In particular, substituting n = in (jSJ) gives ([!]) and putting n — 1 yields (JS} . It is easily 
seen that for n > 1, formula (jHJ) contains fewer summations than the corresponding 
formula for ((4n + 3) given by (j3J). 

There exists a bivariate unifying formula for identities ([2]) and (j7|) 

_l)fc+i 5P-x 2 ^(m 2 -x 2 ) 2 + V 



f fc If (-l) w 5^-^ W 

KJ ^k A -x 2 k 2 -y 4 2^ k( 2k ) k A - x 2 k 2 -y 4 11 m 4 - x 2 m 2 -y A ~ 

fc=l y fc=l Mfc/ y m=l y 

It was originally conjectured by H. Cohen and then proved by D. Bradley p] and, 
independently, by T. Rivoal [22]. Their proof consists of reduction of ([9]) to a finite 
non-trivial combinatorial identity which can be proved on the basis of Almkvist and 
Granville's work [1]. Another proof of ([9]) based on application of WZ pairs was given 
by the authors in [12]. Since 



< 10 > E fc4 _ _ „4 = E E (" + m )c(2r» + 4m + 3). V", |»|» + |v|*<l, 



n=0 m=0 



the formula ([9]) generates Apery-like series for all ((2n + 4m + 3), n, m > 0, convergent 
at the geometric rate with ratio 1/4 and contains, as particular cases, both identities 
@ and (J7|). Indeed, setting x = a and y = yields Koecher's identity (T5]), and setting 
x = 0, y = a yields the Borwein-Bradley identity (JTj). Putting 



2 x 2 + v/^TV , 2 x 2 -^TW 

11 a := , o := . 

2 ' 2 



we can rewrite (Q in a more symmetrical way 

, , f fc 1 f ( Z 1)^W - a 2 - 6 2 )(1 ± a ± &)„_! 

1 j 2-(P-a 2 )(P-6 2 ) 2 2_. n( 2 ;)(l±a)„(l±6) n 

Here and below (w ± v ± tu) means that the product contains the factors u + v + w, 
u + v — w,u — v + w,u — v — w. 

In [12], the authors showed that the generating function fflOl) also has a much more 
rapidly convergent representation, namely 

(u) f k _ i f (-ir-v(n) nr=\(K - * 2 ) 2 + V) 

where 

r(n) = 205n 6 - 160ri 5 + (32 - 62x 2 )n 4 + 40xV + (x A - 8x 2 - 25y 4 )n 2 + 10y 4 n + y\x 2 - 2). 
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The identity ( TT3"j) produces accelerated series for all (,{2n + 4m + 3), n, m >, convergent 
at the geometric rate with ratio 2~ 10 . In particular, if x = y = we get Amdeberhan- 
Zeilberger's series [2] for C(3), 

(M) c(3) = 1 (- 1 ) n ~ 1 (^5n 2 - 160n + 32) 



2^ n 5M 5 



n=l 



It is worth pointing out that both identities (J9J) and ([13]) were proved in [12] by using 
the same Markov-WZ pair (see also [TJ1 p. 702] for the explicit expression), but with 
the help of different summation formulas. 

A more general form of the bivariate identity for the generating function 

OO OO / \ 

S ( m J ( A o^ 2n + 4m + 4) + 5 C(2n + 4m + 3) + C C(2n + 4m + 2))x 2n y im 



n=0 m=0 

.2 



2^ J,4_ x 2 k 2_ 4 ' l X l + M <X ' 

fe=l y 

where Ao, -Bo, Co are arbitrary complex numbers, was proved in [12] by means of the 
Markov- Wilf-Zelberger theory. More precisely, we have 

v^ A + .Bpfe + CpA: 2 = y> dn 

V 1 ^ i-4 _ ^,2i,2 _ „4 TT™ f 



^-^ A; 4 — x 2 &; 2 — y 4 ■ / — J' rim=i( m4 — x 2 m 2 — y 4 



where 



n-l 



(-lf-^O^n 2 - X 2 ) yj 2 2 2 4 , 

In/ m=l 

(40n + !0)L n + (35n 5 - 35n 3 x 2 + 4n(3a; 4 + 10?/ 4 ))L n _i 
+ 4(5n 2 - 2x 2 ) 

and L n is a solution of a certain second order linear difference equation with polynomial 
coefficients in n and x, y with the initial values L Q = C , L\ = (5 — 2x 2 )A /15 + (5x 2 — 
1 - 4(x 4 + 6y 4 ))C /30. If we take A = C = 0, B = 1 in QX5J, then L n = for all 
n > and we get the bivariate identity (JS}. 

First results related to generating function identities for even zeta values belong to 
Leshchiner [16] who proved (in an expanded form) that for \a\ < 1, 
(16) 

_ 1 \ r( 2n_^ _lf 1 3A: 2 + a 2 / _ a 2 ^ 

^-V 1 2 2 ™+ 1 J C(2n + 2)a n 2 -a 2 ~ 2^^( 2fc ) k 2 - a 2 H l* m 2 J ' 

n=0 v 7 n=l fe=l ^ V fc / m=l v 7 

Comparing constant terms on both sides of (1161) yields 



In 2006, D. Bailey, J. Borwein and D. Bradley [5] proved another identity 



(17) Y((2n + 2)a 2n = Y — i— = 3V TT ( 



fc— 1 / , ■) 

m — 4a 



r(t)(^ 2 -« 2 )^ A A- 



2 - a 2 



It generates similar Apery-like series for the numbers ((2n + 2), which are not covered 
by Leshchiner's result ( IT6|) . In the same paper [5], a generating function producing 
fast convergent series for the sequence ((2n + 4), n = 0,1,2, ... , was found, which for 
\a\ < 1, has the form 



here B 2n G Q are the even indexed Bernoulli numbers generated by 

xcoth(x) = ^5 2n — — . 

n=0 ^ ' ' 



i " 2 

2 



It was shown that the left-hand side of (118]) represents a Maclaurin expansion of the 
function 

na csc(na) + 3 cos(7ra/3) — 4 
4a 1 ' 
Comparing constant terms in ( fl8|) implies that 

36^ 1 



A': 



The identity f fT8|) gives a formula for £(2n+4) which for n > involves fewer summations 
then the corresponding formula generated by ffl~6j) . Note that a unifying formula for 
identities generating even zeta values similar to the bivariate formula OH]) for the odd 
cases is not known. 

The Hurwitz zeta function defined by 



C(«,*) = £ 



k=0 (k + vy 

for s G C, Res > 1 and v ^ 0,-1,-2,... is a generalization of the Riemann zeta 
function ((s) = ((s, 1). In this paper, we prove a new identity for values ((2n, v) which 
contains as particular cases Koecher's identity (J2]), the Bailey-Borwein-Bradley identity 
(TTT1) . some special case of identity (|9]) and many other interesting formulas related to 
values of the Hurwitz zeta function. We also get extensions of identities ([9]) and ( TT3|) 
to values of the Hurwitz zeta function. The main tool we use here is a construction of 
new Markov- WZ pairs. As application of our results, we prove several conjectures on 
supercongruences proposed by J. Guillera and W. Zudilin pi)], and Z.-W. Sun [251 126] . 
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2. Background 



We start by recalling several definitions and known facts related to the Markov- Wilf- 
Zeilberger theory (see [T7J [HJ H]). A function H(n, k), in the integer variables n and 
k, is called hypergeometric or closed form ( CF) if the quotients 

H(n + l,k) , H(n,k + 1) 

and 



if(n,fc) H(n,k) 

are both rational functions of n and k. A hypergeometric function that can be written 
as a ratio of products of factorials is called pure-hypergeometric. A pair of CF functions 
F(n, k) and G(n, k) is called a WZ pair if 

(19) F(n + 1, k) - F(n, k) = G(n, k + 1) - G(n, k). 

A P-recursive function is a function that satisfies a linear recurrence relation with 
polynomial coefficients. If for a given hypergeometric function H(n,k), there exists a 
polynomial P(n, k) in k of the form 

P(n, k) = ao(n) + ai(n)k + ■ ■ ■ + aL(n)k L , 

for some non- negative integer L, and P-recursive functions ao(n), . . . , a^(n) such that 

F(n, k) := H(n,k)P(n,k) 

satisfies f fT9|) with some function G, then a pair (F, G) is called a Markov- WZ pair 
associated with the kernel H(n,k) (MWZ pair for short). We call G(n,k) an MWZ 
mate of F(n, k). If L = 0, then (F, G) is simply a WZ pair. 

In 2005, M. Mohammed [18] showed that for any pure-hypergeometric kernel H(n, k), 
there exists a non-negative integer L and a polynomial P(n, fc) as above such that 
F(n, k) = H(n, k)P(n, k) has an MWZ mate G(n, k) = F(n, k)Q(n, k), where Q(n, k) 
is a ratio of two P-recursive functions. Paper [19J is accompanied by the Maple package 
MarkovWZ which, for a given H(n, k) outputs the polynomial P(n, k) and the G(n, k) 
as above. 

From relation (Tl9|) we get the following summation formulas. 
Proposition A. [TBI Theorem 2(b)] Let (F, G) be an MWZ pair. If lim F(n, k) = 

for every k > 0, then 

oo oo oo 

(20) V F(0, k) - lim ]T G(n, fc) = V G(n, 0), 

i - — ' fc— s>oo ' ^ — ' 
fc=0 n=0 n=0 

whenever both sides converge. 

oo 

Proposition B. [TS1 Cor. 2] Lei (F, G) be an MWZ pair. If lim £ G(n, fe) = 0, t/ien 

fe-^oo n=0 

oo oo 

(21) £ F(0, fe) = ]T(F(n, n) + G(n, n + 1)), 

fe=0 n=0 

whenever both sides converge. 

Formulas ( 120]) , (}2~Tj) with an appropriate choice of MWZ pairs can be used to convert 
a given hypergeometric series into a different rapidly converging one. 



To ensure wider applications of WZ pairs for proving hypergeometric identities we 
use an approach due to I. Gessel [5] (see also [201 §7.3, 7.4]). It is based on the fact that 
if we have a WZ pair (F, G), then we can easily find other WZ pairs by the following 
rules. 

Proposition C. [8, Th. 3.1] Let (F, G) be a WZ pair. 

(i) For any complex numbers a and (3, (F(n + a, k + (3), G{n + a,k + (3)) is a WZ 
pair. 

(ii) For any complex number^, ( , fF(n,k),'yG(n,k)) is a WZ pair. 

(iii) If p(n, k) is a gamma product such that p(n + 1, k) = p(n, k + 1) = p(n, k) for 
all n and k for which p(n,k) is defined, then (p(n,k)F(n,k),p(n,k)G(n,k)) is a WZ 
pair. 

(iv) (F(—n,k),—G(—n — l,k)) is a WZ pair. 

(v) (F(n, -k), -G(n, -k + 1)) is a WZ pair. 

(vi) (G(k,n), F(k,n)) is a WZ pair. 

The WZ pairs obtained from (F, G) by any combination of (i)-(v) are called the 
associates of (F, G). The WZ pair of the form (vi) and all its associates are called the 
duals of (F, G). 

3. The identities 



Theorem 1. Let a, a, (3 G C, \a\ < 1, a ^ (3, and a ± a, ± a be distinct from 
0, —1, —2, .... Then we have 



a 2n (C(2n + 2, a) - ((2n + 2, 0)) 



(3 — a 



n=0 



(22) = y (-l)"-i(l + g - (3)^(1 + (3- q)n_i(l + 2a) n _ 1 (l - 2a) n _! 

(a + a) n (a - a) n (/3 + a) n (/3 - a) n 

(5n 2 + 3n(a + (3-2) + 2(a- l)((3 - 1) - 2a 2 ) 

x <T) ' 

Proof. By the definition of the Hurwitz zeta function, we have 

Now define a Markov kernel H(n, k) by the formula 

jj, js _ (at + a)k{a - a)k((3 + a)k((3 - a)k(n + 2k + a + f3) 
(a + a) n+k+ i(a - a) n+k+ i((3 + a) n+k+ i((3 - a) n+k+ i 
Applying the Maple package Markov-WZ we get the associated WZ pair 

(_1)«(1 + a _ /3 ) n (i + p _ a ) n (i + 2o ) n (i - 2a) n 



F(n, k) = H(n,k) 
G(n,k) = F(n,k) 



( 2 D 

5n 2 + n(3a + 3/3 + 4) + 2a(3 - 2a 2 + (2k + 1)(1 + a + (3) + 2k(k + 3n) 

2(2n+ l)(n + 2k + a + (3) ' 
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Now by Proposition A, we obtain 

oo oo 

^F(0,fc) = ^G(n,0), 

k=0 n=0 

which implies fl22l) . □ 

Multiplying both sides of (1221) by f3 — a and letting (3 tend to infinity we get an 
extension of the Bailey- Borwein- Bradley identity to values of the Hurwitz zeta function: 

(24) ± a» X (2n + 2, a) = £ < 3 " + 2 ° ~ 2 > (1 + ^ F ^ • 

Setting a = 1 in (124]) yields the Bailey-Borwein-Bradley identity ( fT7|) . 

Replacing a by a/2 and a, (3 by 1 + a/2, 1 — a/2, respectively, in (122]) and taking into 
account (|23l) . we get Koecher's identity (T5]). 

Letting /3 tend to a in (122]) and using the equality 

we get the following. 

Corollary 1. Let a, a e C, |a| < 1, ana 1 « ± a ^ 0, —1, —2, .... Then 

oo OO , 

g(n + l) C (2n + 3,a)^ = g ((j , +a ; 2 a _ a2)2 

k 

1 (-ly^jhn 2 + Qn(a - 1) + 2(a - l) 2 - 2a 2 ) (n - 1)! 2 (1 + 2a) w _i(l - 2a) w _! 
"2^ ^) (« + «) 2 (a-a) 2 " 

Taking a = 1 in Corollary [T] we get the following identity for odd zeta values. 

Corollary 2. Let a G C, \a\ < I. Then 
(25) 

f nC(2„ + l)a— - f * - I f (-l)-'(5n'-<.') fi 1-4,'M' 



Note that the right-hand side equality of (125]) also follows from the bivariate identity 
(Q or (fT2l as was shown by D. Bradley (see [7] Cor. 1]). 

It is clear that identity f]25]) gives formulas for odd zeta values which are linear com- 
binations of series generated by the bivariate identity (JHJ). Thus comparing constant 
terms on both sides of (125 p gives Apery's series ([T]) for C(3)- Similarly, comparing coeffi- 
cients of a 2 gives formula (jl]) for ((5). It produces the following complicated expression 
forC(7): 

u " 6 h k7 i 2 t) *k H 2 t) k? 2 Q k 



10 " (-i)fc-i ^i^i 
+ 3 ^ k 3 ( 2k ) 2s f 2s m 2' 

k=l \kJ j=l J m=l 
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which can be written as 

C(7) = i(4K--B), 

where K and B are right-hand sides of formulas (jSJ) and ([6]), respectively. More gener- 
ally, if we denote 

fc-i 

gU(k) :=[t r ]l[(l-r s t)- 2 , 

3=1 

then taking into account that 

5fc 2 - 2a 2 _ 2 3 _ , 

(A; 2 — a 2 ) 2 ~ l-a 2 /k 2 + (l-a 2 /k 2 ) 2 ~ f^ Q k 2 i 

and comparing the coefficients of a 2n on both sides of ( 1251) . we get 
Corollary 3. Le£ n be a non-negative integer. Then 

(26) C(2n + 3) = -— J> + 5) £ £(-4)^*) 

"'" j=0 fc=l UJ r=0 

Consider several other particular cases of Theorem HJ Replacing a by a/2, a by 1/2, 
and /3 by 1 in (122]) and noting that 

y, 2fc + 3/2 = y (-l)"' 1 

^ ((jfe + 1/2) 2 - a 2 /4)((£; + l) 2 - a 2 /4) ~ ^ n 2 - a 2 ' 
we get the following identity. 

Corollary 4. Lei a be a complex number, distinct from a non-zero integer. Then 
;_l)n-i ! ~ (_i)n-i(io ra 2 -3n-a 2 ) 



4^ fl (2«-l)(^-^)( 2 ;)n; =1 (l-aV(n+j) 2 )' 
In particular, 

C(2) - 1 V (-^""'(IQ" ~ 3) 

2 £t - 2 (2--i)( 2 ;) ' 

Substituting a = l/3, /3 = 2/3, a = 0in Theorem [I], we get 
Now observing that 

(2\ _ (3n)! 

and 



n \ / n 



3 7 „ V 3 / „ 27 n n! 



CO 



((2,l/3)-((2,2/3) = 9^^=: 9£T 



n=l 



<) 



(where (^) is the Legendre symbol), we get the following formula 



n=l 



(15n-4)(-27) 

n3 {n) („) 



ra-1 



which was conjectured by Z.-W. Sun in [25] . 

Substituting a = 1/4, — 3/4 recovers Theorem 3 from [T3] and in particular (when 
a = 0), it gives the following formula for Catalan's constant G := Xl^Lo (2n+i) 2 : 



1 ^> (-l) n - 1 256 n (40n 2 -24n + 3) 



G=-Y 

fizL Z-w 



64^ (4n\2(2n 



n=l 



o ( 2 ;)- 3 (2--i) 



Applying Proposition B to the Markov-WZ pair found in the proof of Theorem [Tj, we 
get the following identity which generates Apery-like series for the differences ((2n + 
2, a) — ((2n + 2, 0) converging exponentially fast as 2~ 10 . 

Theorem 2. Let a, a, (3 G C, \a\ < 1, a ^ 0, and a ± a, ± a be distinct from 
0,-1,-2 Then 



a — J2^ 2n + 2 . «) - C(2^ + 2, /3))a 2n 
— a 



(27) n= ° 
y> (-1)"- Wn) ~ (<» ~ OT(j 2 - 4a 2 ) 

where 

PaA n ) = 2(2^ -l){3n + a + 0- 3)((2n - 1 + «) 2 - a 2 )((2n - 1 + /3) 2 - a 2 ) 
+ ((n + a - l) 2 - a 2 )((n + - l) 2 - a 2 )(13n 2 + 5n(a +0-2) + 2((1 - a)(l - /3) - a 2 )). 

Setting a = 1/3, /3 = 2/3, a = in Theorem [21 we get the following fast converging 
series for the constant K : 

_ (-27) n_1 (5535n 3 - 4689n 2 + lllOn - 80) 

n3(3n-l)(3n-2)e) 2 ( 3 n ") ' 

Setting a = 1/4, = 3/4, we recover Theorem 4 from [T3] . 

Multiplying both sides of (12~7]) by /3 — a and letting /3 tend to infinity we obtain an 
extension of Theorem 2 from [TTJ to values of the Hurwitz zeta function 

(28) f C(2n + 2, a)a^ = V p %\ ^M^lI , 

where p(n) = 2(2n - l){{2n - 1 + af - a 2 ) + (5n + 2a - 2)((n + a - l) 2 - a 2 ). In 
particular, setting a = 1, a = in (128]) we get Zeilberger's series [301 §12] for C(2), 

21n-8 



c(2) = E- 

n=l n 



Replacing a by a/2 and a, by 1 + a/2, 1 — a/2, respectively, we recover Theorem 4 
from [TT]. Letting tend to a in (1271) . we get the following 
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Corollary 5. Let a, a e C, \a\ < 1, and a ± a fre distinct from 0, — 1, —2, .... Taen 

oo 



)a 2k - 2 



k=l 



1^ l-4a 2 /m 2 



__/2n\ 5 //_ , _ i\Q _9\9//„ i _ 1\5 _9. I a\-> AA 



2 ^ n( 2 n ") 5 ((n + a - If - a 2 ) 2 ((n + a - l) 2 - a 2 /4) 2 JLi ((1 + ^) 2 - (^) 2 
u>aere 

p a (n) := p a ,a( n ) = 2 ( 2n - ^C 3 " + 2a - 3)((2n + a - l) 2 - a 2 ) 2 

+ ((n + or - l) 2 - a 2 ) 2 (13n 2 + 10n(a - 1) + 2((1 - a) 2 - a 2 )). 

Setting a = 1 in Corollary we get the following identity. 

Corollary 6. Let a e C, |a| < 1. Taen 

h 2^ n(2:) V-a^-aV4) 2 i 1 1 (l-« 2 /(- + -) 2 ) 2 ' 

where p(n) = 2(2n - l)(3n - l)(4n 2 - a 2 ) 2 + (n 2 - a 2 ) 2 (13n 2 - 2a 2 ). 

Setting a = in Corollary [6] we get Amdeberhan-Zeilberger's series ( I14p for C( 3 )- 
The next theorem gives a generalization of identity f lT2|) . 

Theorem 3. Let a,a,b G C and a ± a, a ± b be distinct from 0, —1, —2, .... Taen the 
following identity holds: 

oo , 

Ere + a 
((A; + a) 2 -a 2 )((A; + a) 2 -& 2 ) 

(29) k ~° 

1 (-l)"'^! ±a± 6)n-i(5n 2 - 6n(l - a) + 2(1 - a) 2 - a 2 - 6 2 ) 

2 ~f n( 2 n n ) (a ± a)„(a ± 6) n 
Proof. Taking the kernel 

(a + a)fc(a — a)fc(a + o) fe (a — o)&(n + 2/c + 2a) 



H(n, k) 



F ( n > k ) = S^H 1 ±a± b )nH(n, k) 



(a + a) n+k+ i(a - a) n+k+1 (a + b) n+k +i(a - b) n+k+1 
and applying the Maple package MarkovWZ we get that 

-iy 

o 

and 

IN „, , x 5n 2 + 6an + 4n + 2a + 2a 2 + l-a 2 -o 2 + H2fc + 6n + 4a + 2) 

Gin, k) = Fln.k) ; ; : 

V ' ; V ' ; 2(2n + l)(n + 2A; + 2a) 

give a WZ pair, i.e., 

F(n + 1, k) - F(n, k) = G(n, fc + 1) - G(n, k). 

n 



Now by Proposition A, we get 

oo oo 

£>((),*) = £G(n,0), 

k=0 n=0 

which implies fl29l) . □ 

Making the substitution (TTTT) in (I2"§j) we get a generalization of Cohen's identity to 
values of the Hurwitz zeta function. 

Corollary 7. Let x,y,a G C, \x\ 2 + |y| 4 < 1 and —1, —2, Then 

OO j OO OO / \ 

£ (^-^+^ ■ ££ (" - ") c(2n + 4m + 3 

1 ^ (-l) w ~ 1 (5n 2 - 6w(l - a) + 2(1 - a) 2 - x 2 ) l4 (j 2 - x 2 ) 2 + % 4 



2^ n( 2 n ") ((n + a - l) 4 - x 2 (n + a - I) 2 - y 4 ) f = \ {j + a- l) 4 - x 2 (j + a - l) 2 - y 4 
Setting a — 1/2, x = y = in Corollary [7] we get the following formula: 



1 ^ (-l)"- 1 (10n 2 -6n + l)256 n 

Y8 k ^Sf ' 



Applying Proposition B to the Markov- WZ pair used in the proof of Theorem [3] we get 
the following identity. 

Theorem 4. Let a, a, b G C and a ± a, a ± b ^ 0, —1, —2, .... TTien 

+ a 



fc=0 



((& + a) 2 -a 2 )((£; + a) 2 -6 2 ) 



1 (-1)" X (l ± a ± 6)n-l(Qi ± fl)n-l(Q! ± &)n-l , s 

2 h n( 2 :)(a±a) 2n (a±b) 2n 
where 

q(n) = 2(2n - l)(3n + 2a - 3)((2n + a - l) 2 - a 2 )((2n + a - l) 2 - 6 2 ) 

+ ((n + a - l) 2 - a 2 )((n + a - l) 2 - 6 2 )(13n 2 - 10n(l - a) + 2(1 - a) 2 - a 2 - 6 2 ). 

Making the change of variables (ITTj) in Theorem H] we get a generalization of the 
identity ( fT3l) to values of the Hurwitz zeta function. 

Corollary 8. Let x,y,a G C, \x\ 2 + \y\ 4 < 1 and a —1, —2, .... Then 



(30) 



iA(-i)-%) n;:i((j 2 -^ 2 ) 2 +% 4 ) 

2 ^ n; =0 ((^ + j + « - 1) 4 - ^ + j + « - i) 2 - ?/ 4 ) 
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where 

Q{n) = 2(2n - l)(3n + 2a - 3)((2n + a - l) 4 - x 2 (2n + a - l) 2 - y 4 ) 

+ {(n + a- l) 4 - x 2 (n + a - I) 2 - y 4 )(13n 2 - 10n(l - a) + 2(1 - a) 2 - x 2 ). 

Formula ( 1301) produces accelerated series for the values C(2^ + 4m + 3, a), n,m > 0, 
convergent at the geometric rate with ratio 2~ 10 . 

Remark. Note that Theorem H] can also be obtained from Proposition A applied to 
the Markov- WZ pair associated with the kernel 

f oi\ tt i U \ (a + a) k+n (a - a) k+n (a + b) k+n (a - b) k+n , ^.o ^ 

(31) H 1 (n,k) = - r r (3n + 2k + 2a). 

{a + a) 2 n+ k +i{a - a) 2n+k+ i(a + 0) 2n+k+ i{a - o) 2n+k+ i 



4. Hypergeometric reformulations 



Theorems (TJ [3] admit nice reformulations in terms of the hypergeometric and digamma 
functions which can be useful in their own. Let us recall the definition of the digamma 
function that is the logarithmic derivative of the Gamma function 

d °° 

*l){z) = — iogr(z) = -i + J2 



dz 



n=0 



n 



n + z 



Obviously, = —7, where 7 is Euler's constant. The function ip(z) is single-valued 
and analytic in the whole complex plane except for the points z = —m, m = 0, 1, 2, . . . , 
where it has simple poles. It's connection with the Hurwitz zeta function is given by 
the formula 

d k 

^ k \z) = —^(z) = (-l) k+1 kK(k+l,z). 

Now it is easily seen that Theorems [1] and [3] are equivalent to the following statements. 

Theorem 5. Let a, a, (3 be complex numbers such that a ± a, (3 ± a be distinct from 
0, —1, —2, .... Then we have 

a(a - /3)((1 + a)(l + (3) + a(3 - 2a 2 ) 



ip((3 + a) — ip(a + a) + ip(a — a) — i/j(/3 — a) 



a 2 W 



x S F 7 



1, 1, 1 + a - p, 1 + - a, 1 ± 2a, | + ± ^ 



§, 1 + a ± a, 1+ P±a, § + 



10 

3(a+/3) I VD 

10 10 



10 



where D = 40a 2 + 9(a - P) 2 - 4(a - - 1). 

Theorem 6. Le£ a, a, 6 fre complex numbers such that a ± a, a ± 6 6e distinct from 
0, —1, —2 T/ien we /iawe 



ij){a + b) + i[)(a — b) — ip{a + a) — ip(a — a) 
/ 



(a 2 -6 2 )((l + «) 2 + « 2 



x 8 F 7 



l,l,l±a±6,^f2± 



2(a 2 -a 2 )(a 2 -o 2 ) 



¥) 2 



1 + a ± a, 1 + a ± 6, 
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2+3« _(_ 
5 nZ 



a?+b 2 
5 



/ 



Dividing both sides of the identity in Theorem by (3 — a and letting (3 tend to a, 
we get the following. 

Corollary 9. Let a, a be complex numbers such that a±a be distinct from 0, —1, —2, .... 
Then we have 



ijj'ia — a) — ijj'ia + a) 



a((l + a) 2 + a 2 - 2a 2 
(a 2 - a 2 ) 2 



x «F 7 



l,l,l,l,l±2a,^± v / 2 fi-(i^)2 



V 



|, 1 + a ± a, 1 + a ± a, ^ ± - (±f^) 2 



Multiplying both sides of the above equality by (a — a) 2 and letting a tend to a, we 
get the following identity. 

Corollary 10. Let a be a complex number such that 2a be distinct from 0, —1, —2, .... 
Then we have 



1,1,1 -2a. 



7+3a±V9a 2 +2a-l 



5^4 



2 - 



l + 2a, 



4 



4a 



l + 2a 



Taking, for example, a = 1/4 in the above identity we get the following non-trivial 
summation formula: 

^ (-l) fc - 1 (5fc-2) L 



k=l 



Corollary 11. Let n be a non-negative integer. Suppose that a, (3 are complex numbers 
such that 2f3 and a + (3 are distinct from —n, —n — 1, — n — 2, ... , and a — (3 ^ n,n — 
1, n — 2 Then we have 



7-^6 



(32) 



n + 1, n + 1 ± (a - fi), n + 1 ± (2n + 2/3), n + | + 
§ + n, n + a + 1 ± (/3 + n), 2n + 2/3 + 1, n + f + 



7 , 3(a+/3) 



10 



± 



3(a+/9) I 

10 10 



/Pi 
10 



4 



(n + a + /3) n +i(« + 2) 

n+l 



(l + a-/3) n+1 (l + 2n + 2/3) n+1 ' 
w/iere D a = 40((n + /3) 2 - (a - l)(/3 - 1)) + 9(a + /3 - 2) 2 . 

Proof. To prove ( 132|) . it is sufficient to consider both sides of the identity from Theorem[5] 
as meromorphic functions of variable a and compare corresponding residues at the 
simple pole a = (3 + n, where n is a non-negative integer, on both sides. Imposing 
the restrictions on a, (3 formulated in the hypothesis implies that the residue on the 
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left-hand side of the identity of Theorem is equal to —1, and clearly, we have 

a(a - /3)((1 + a)(l + /3) + a/3 - 2a 2 ) 



1 = lim (n + j3 — a) 



la* 



a?)(ff- 



E 



fc!(l ±{a- P)) k {l ± 2a) fc (7/5 + 3(a + /3)/10 ± \AD/10) 



, n (-4) fc (3/2) fc (l + a ± a) fc (l + /3 ± a) fc (2/5 + 3(a + /3)/10 ± v^/10)* 
(n + /3)(a - /?)((! + a)(l + /3) + a/3 - 2(n + /3) 2 ) 



E 



(a 2 - (n + /3) 2 )(n + 2/3)(-l)™ 

fc\ (1 ± (a - /3)) fc (l ± (2n + 2/3))*((14 + 3(a + /3) ± v^/lO)* 

n ' 



. , 7 x -4) k (3/2) k (l + a±(/3 + n)) k (l + 2/3 + n) k ((4 + 3(a + (3)±VTT 1 )/10) k 
(n + /3)(a - /3)(1 ± (a - /3)) n (l ± (2n + 2/3)) n (10n + 4 + 3(a + /3) ± 



20(a 2 - (n + /3) 2 )(n + 2/3)4"(3/2) n (l + a ± (/3 + n)) n (l + 2/3 + n), 



x 7 F 6 



'n + 1, n + 1 ± (a - /3), n + 1 ± (2n + 2/3), n + J + ± 
| + n, 7i + a + 1 ± (/3 + n), 2n + 2/3 + 1, n + § + ± ■ 



10 
/237 



10 



4 



Now taking into account that 
a-/3 (1 + P-a) 



-iy 



a — (3 — n (1 + a — (3 — n) n 

(10n + 4 + 3(a + /3)) 2 - D 1 
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n + P (1 - 2n - 2/3) n _ (-1)' 
n + 2/3 ' (l + 2/3 + n) n ~ 2 

(a-/3 + n + l)(3n + 2/3 + l), 



□ 



and 4 n (3/2) n = (2n + l)!/n!, we get the required identity 
Similarly, from Theorem [6] we get 

Corollary 12. Let n be a non-negative integer and a, a be complex numbers such that 
2a 7^ —n, — n — 1, ... , and a ± a ^ —n — 1, — n — 2, ... , — 2n — 1. Then we have 



( 



(33) 



n + l,2n + a + l±a, l-a±a, n + ± 



a 2 +(n+a) 2 _ ^ l-a ^2 



+ n, 2n + 2a + 1, n + a + 1 ± a, n + ^ ± J ' *±Ss±°i£ - (1=2)2 



(n + 2) n+1 (n + 2a) 



ra+1 



(n + a + 1 ± a) n+ i 
Proof. To prove (133]) . we rewrite the identity of Theorem [6] in the form 



E 



k + a 



((1 + a) 2 + a 2 - a 2 - b 2 ) 



(34) 



, o ((k + a) 2 - a 2 )((k + a) 2 - b 2 ) ~ A(a 2 - a 2 )(a 2 - b 2 ) 
1, 1, 1 ± a ± b, ± v/1^6 2 ) 



|, 1 + a ± a, 1 + a ± b, 2±ps 
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4 



where 



D{¥) : 



2 , a 2 + Ir 



1 — a 



replace b 2 by z and consider both sides of (I34p as meromorphic functions of variable z. 
Suppose that 2a is distinct from — n, —n — 1, .... This restriction ensures that (n+a) 2 ^ 
(J + a) 2 for any non-negative integer j ^ n. Now equating residues on both sides of 
(I34p at the simple pole z = (n + a) 2 , where n is a non- negative integer, we have 



n + a 



(n + a) 2 



lim ((n + a) 2 — z) 



z— >(n+a) 2 



(1 + a) 2 + a 2 - a 2 - 
A(a 2 - a 2 )(a 2 - z) 



x 



E 

k=n 



n 



(0'±«) 



(-4)*(3/2) fc (^ ± /D(^)) fc A = l ((j + a )2 _ z )^j + a )2 _ a 2) 

(1 ± a ± (n + a))k 



(1+Qf " a2 "' ! " 2na E / '' 



(-l)»-i(2a + n) n 



k=n 



-4) fc +!(3/2) fe (a ± a) fc+1 (2a + 2n + l) fc . 



± v/D((n + a) 2 )) fe _ (1 ± a ± (n + a))„((2n + a + l) 2 - a 2 



(2±pa ± v /£,(( n + a )2)) fc 4™+ 1 (3/2)„(2« + n) n (a ± a) n+1 



x 7-Ffi 



n 



+ 1, 2n + a + 1 ± a, 1 - a ± a, n + ± ^((n + a) 2 ) 



2n + 2a + 1, n + a + 1 ± a, n + ^ ± y/D((n + 
Here in the last equality we used the fact that 

± y/D((n + a) 2 )) n (2n + a + l) 2 - a 2 



a) 



4 



' 7+3a 
, 5 



( 2 ^±V^((^ + «) 2 ))„ (l + «) 2 -a 2 -n 2 -2n«' 

Finally, after simplifying based on the cancelations 

(l-a-n-a) n _ (-1)" (1 + a - n - a) n _ (-1)" 

(a + a) n+ i a + a + n ' (a — a) n +i a — a + n ' 

we get the required identity. 

Setting a = 1 and replacing n by n — 1 in (|33|) we get 



□ 



Corollary 13. Lei n be a positive integer and a be a complex number such that ±a ^ 
—n — 1, — n — 2, ... , — 2n. Then we have 



( 



(35) 



7-^6 



n, 2n ± a, ±a, n + 1 ± 



a 2 +n 2 



n + |, 2n + 1, n + 1 ± a, n ± a/ 



/ 



2n 

n 



j 2 — a 2 



Formula (135]) generalizes similar identities from [7J §2]. In particular, substituting 
a = \fcjn gives [TJ Corollary 2] and substituting a = iy/b + n 2 gives Corollary 3]. 
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5. SUPERCONGRUENCES ARISING FROM THE AMDEBERHAN-ZEILBERGER SERIES 

FOR C(3) 

In this section, we consider supercongruences arising from the Amdeberhan-Zeilberger 
series f jl4"|) for C(3). In [ID] . Guillera and Zudilin proposed conjecturally ap-adic analog: 

^y^|(205fc 2 + 160fc + 32)(-l) fc 2 10fc = 32p 2 (mod p 5 ) for prime p > 3. 

In [26], Z.-W. Sun formulated more general conjectures: let p be an odd prime, then 
p-i /2k\ 5 

^(205fc 2 + 160fc + 32)(-l) fe ( J = 32p 2 + 64p 3 F p _i (mod p 7 ) forp ^ 5, 

where = £)fc=i V^j anc ^ 

(p- 1 )/ 2 /2Jfc\ 5 896 

(36) ^ (205/c 2 + 160A;+32)(-l)M J = 32p 2 + — p 5 5 p _ 3 (mod p 6 ) forp > 3, 

A;=0 ^ ' 



where Bq, B\, B2, . . . are Bernoulli numbers. 

Moreover, Sun [26] introduced the related sequence 



I n ~ 1 /9k\ 1 



fc=0 



and conjectured that for any positive integer n, a n should be a positive integer. 

In this section, we confirm these conjectures (with the only exception that we prove 
modulo p 5 ) and prove the following theorems. 



Theorem 7. Let n be a positive integer and let 

n-l /„, x 5 



A n := ^(-l)"- 1 ^ ( ' J (205A; 2 + 160A; + 32). 



k=0 

Then the following two alternative representations are valid: 

n-l 



k=0 

and 

2 "--1 /r> 1 \ /r> ; r>\ 2 



Theorem [7] implies immediately the following 



Corollary 14. For any positive integer positive integer. 
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Theorem 8. Let p be an odd prime. Then the following supercongruences take place: 
^(205fc 2 + 160A; + 32)(-l) fc ( J = 32p 2 + 64p 3 H p _ 1 (mod p 7 ) forp^ 5, 



k=0 

(P-l)/2 /m\5 



/2k\ 

(205A; 2 + 160A; + 32)(-l) /c ( J = 32p 2 (mod p 5 ) for p > 3. 

k=0 ^ ' 

The proof of Theorem [7] is contained in the following two lemmas. 
Lemma 1. For any positive integer N, the following identity holds: 

Proof. Consider the Markov kernel H\{n, k) defined in ( 13 ip for a = l,a = b = 0: 

Then the corresponding to it WZ pair has the form 

(-l)"n! 6 (-l)"n! 6 (fc + n)! 4 

F(n, fc) = — ; — n — HAn, k), Gin, k) = — — — rrrQin, k), 

K ' ; (2n)! 2(2n + 1)!(A; + 2n + 2)! 4 ^ v ' ; ' 

where 

q(n, k) = {n + l) 4 (205n 2 + 250n + 77) + A;(254 + 344A; + 1526n + 3628n 2 + 888k 2 n 

+ 2928kn 2 + 4268n 3 + 248A; 2 + lOlfc 3 + 1648kn + 574n 5 + 2486n 4 + 22A; 4 + 2k 5 
+ 664kn 4 + 2288fcn 3 + 408JfeV + 1048JfeV + UlJfeV + 240k 3 n + 26k 4 n). 

It is easy to show (see Proposition C or [201 Ch. 7]) that the pair (F, G) given by 
C( B ,*) = (-l)-t(*)(" + B *- 1 ) 4 (3t + 2n ), 
— , , x */2Jfe\ f(n + 2k- 1)!\ 4 , 

with 

gi(n, k) = 205k 6 + 160k 5 + 32k 4 + 2n 6 + n 5 (4 + 26k) + n 4 (2 + 42k + Ulk 2 ) 

+ n 3 k(16 + 176k + 408k 2 ) + n 2 A; 2 (48 + 368A; + 664A; 2 ) + n£; 3 (64 + 384A; + 574A; 2 ) 
is its dual WZ pair, for which 

(37) F(n + 1, k) - F(n, k) = G(n, k + 1) - G{n, k). 

It turns out for our further proof that it is useful to consider the usual binomial coef- 
ficient (r) in a more general setting, i.e., to allow an arbitrary real number to appear 
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in the upper index of u) , and to allow an arbitrary integer in the lower. We give the 
following formal definition (see [9j §5.1]: 

^^r^ ^ -teger fc > 0; 
k) 1 0, if integer fc < 0. 

After this elaboration we see that G(n, k) is defined for all integers n, k. Rewriting 
F(n, k) in the form 

[^(SC^W), if MO; 
F(n,k) = < 2(n + l) 2 , if fc = 0,n>0; 

[(D, if fc = 0,n<0; 

we can conclude that F(n, k) is well defined for all n, k G Z. Now we can show that 
relation (137)) takes place for all integers n and fc. Indeed, if k < or if k = and 
n < — 1, then all parts in (137)) are zero. If fc = and n > —1, then G(n,k) = 0, 
(5(n, fc + 1) = 2(2n + 3), F(n + 1, k) = 2{n + 2) 2 , F(n, k) = (n + l) 2 , and relation (ETJ) 
is equivalent to the obvious equality 

2{n + 2) 2 - (n + l) 2 = 2{2n + 3). 

If k > and n + k < 0, then F(n, k) = G(n, k) = 0, 

— . (-l) k [2k\ ( n + 2k ^ 4 

F(n + l,fc) = - v ' w 

and 



nn+hk) = ^ J, | a(B + i,*), 



2fc + 2\ /n + 2fc + l x 4 



G(„,t + i) = (-i)'(* + i)^ + 1 J^ B + i + 1 j(» + a n + 3). 

If moreover, n + k + 1 < 0, then F(n + l,k) = G(n, k + 1) = and flUZD holds. If 
n + k + 1 = 0, then G(n, k + 1) = (-l) fe (fc + l) 2 ^ 2 ) and 

F(n + 1, fc) = t^l ^\ fc) = W (2 fc + 2 )(2fc + 1) = G(n, + 1), 



and therefore (137)) holds. If fc > and n + fc > 0, then canceling common factorials on 
both sides of ( 137)) we get the equality 

(n + 2k)%(n + 1, fc) - {n + k + l)V(n, fc) 

= 2(2fc + l)(n + 2fc) 4 (n + 2fc + l) 4 (3fc + 2n + 3) + fc 5 (n + fc + l) 4 (3fc + 2n), 

which can be easily checked by straightforward verification. 

Now let iVeN. Considering relation (137)) at the point (n — N, fc) we have 

(38) F(n + 1- N, fc) - F(n - N, fc) = G(n — N,k + 1) — G(n - iV, fc). 
Summing both sides of (135]) over fc from to iV — 1 we have 

N-l 

(39) ^2(F(n + l-N,k)-F(n-N,k)) = G(n-N,N)-G(n-N,0) = G(n — N,N). 

k=0 
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Now summing over n from to iV — 1 we get 



N-l 



N-l 



J3(F(0, k) - F(-N, k)) = J2 °( n - N > N )- 

k=0 n=0 

Since F(—N, k) = for k = 0, 1, . . . , N — 1, we obtain 

N-l N-l 
k=0 n=0 

or 



T^- l)k \k) ( 205A;2 + 160fc + 32 ) 

k=0 ^ ' 



-1 



iJV-1 



AT 



27V 
N 



N-l 

E 

n=0 



iV + n - 1 



n 



and the lemma is proved. 

Lemma 2. For any positive integer N, the following identity holds: 



(2n + N), 
□ 



Proo/. Let N e Z, N > 0. Put 
(40) 



2 N-l 



fc=0 



2iV - 1\ /2iV - Jfe - 2 



N + k 



N-l-k 



S N :=2^( k + N )\N + 2k + l). 

k=0 ^ ' 



Now rewriting Sjv in the form of a terminating hypergeometric series, we get 



N 



S N = 2(N + 1)J2 



Changing the order of summation and noticing that 



we get 

(41) S N = 2 



2N 
N 



{a)N-k 



(3N + 1) ■ 6 F 5 



-DHa. 



N 



(l-a-N) k 



1, | - *N,-N,-N,-N,-N 
| - |JV, -2iV, -2iV, -2JV, -2JV 



To evaluate the hypergeometric series on the right-hand side of (|4ip . we apply Whipple's 
transformation [29| (7.7)] which transforms a Saalschutian 4^3(1) series into a well- 
poised 7 F 6 (1) series (see [231 P- 61, (2.4.1.1)]): 



4-T3 



(42) 



f — ax — a 2 , dx, d 2 , —N 



x 7 F 6 



_ (g - dx)N(g ~ d 2 ) N 
{9)n{9 - dx- d 2 ) N 

/ - 1, 5/ + 5, 01, a 2 , di, d 2 , -A 7 " 

5/ ' - 5> / ~ a iJ ' - 02, / ~dx,f -d 2 J + N 
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Setting ai = a 2 = —N, d x = -N, d 2 = 1, / = -3iV in (P]), we get 



(43) 



F , 1,1~IN,-N,-N,-N,-N 
6 -| - fiV, -2iV, -2iV, -2JV, -2iV 



(2N + 1) 



(3N + 1)(N+1 
Therefore from (j4"Tl) and fj43|) we get 



■4-^3 



1, -JV, -N, -N 
-2N, -2N, N + 2 



1 . 



/2iVV (2iV + l) 2 / 1, -N,-N,-N 
bN ~ Z{ N] N + l '^ 3 {-2N,-2N,N + 2 



or 



(44) 



5 



2(2iV + l)! 2 (2iV)! 



(iV + l)!iV! 7 



E 



(-jv)i 



fc=0 



-2iV) 2 (iV + 2) fc - 



Replacing Pocchammer's symbols by factorials in ( )44l) we arrive at 



(45) 



^7V = (iV 



+< + \ 2 )^( 



2iV + l\ [2N - k 



N-kJ \ N-k 



Now replacing iV by N — 1 in fj45|) , and using ( 140]) and Lemma [T] we get the required 
identity. □ 

To prove Theorem [HJ we need several results concerning harmonic sums modulo a 
power of prime p. The multiple harmonic sum is defined by 

1 



.finai, a,2, ■ ■ • , a r ; n 



E 



t Q l h" -2 ■ ■ ■ h a r 

l<k 1 <k 2 <--<k r <n v 1 r 



where n > r > 1 and (ai, 02, ... , a r ) G N r . For r = 1 we will also use the notation 

n 1 

:=H{a-n) = Y J T- a and H n := 



k=l 



The values of many harmonic sums modulo a power of prime p are well known. We 
need the following results. 

Lemma 3. p41 Theorem 5.1] Let p be a prime greater than 5. Then 



H p ^ = Hl 3 \ = (modp 2 ), #^ = #^ = (modp). 



(2) _ rr(4) 



Lemma 4. Let p > 5 be a prime. Then 

1 

2 p 



fifi({l} 2 ;p- 1) = -^fifiS (modp 4 ), fifi({l} 3 ;p- 1) = (mod p 2 ) 



fifi({l} 4 ;j9- 1) = (modp). 
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Proof. Since for n > 1 (see, for example, 

H({l} 2 ;n) = 1 -(Hl-H^), 

H({1} 3 ; n) = \{Hl- SH n H^ + 2H®), 

H{{lY- n) = ±(H* - 6H*HW + 8 H n H^ + 3(tff) 2 - 6H^), 
then by Lemma El we get the required congruences. □ 
Lemma 5. Let p > 5 be a prime. Then 



^=1 +pH ^-£hV (modp 5 ). 



Proof. It is readily seen that 

p-i / \ p-i 



l/2p\ = (2p-l 
2\p) \p 



i 1 )=n( 1+ 7)=E^(W;p- 1 ) 



and therefore 

I M = i+ P ff H Vff({i} J ;p-i) Vff({i} 4 ;p-i) (mod P 5 ). 

Now by Lemma HI we get the required congruence. □ 
For similar congruences related to the central binomial coefficients, see [27J. 

Proof of Theorem O From Theorem [7] with n — p, where p > 5 is a prime, we have 
(46) 

For the sum on the right-hand side of f[46p . changing the order of summation we have 

p— 1 /ri 1 \ /r> T r>\ 2 p— 1 



2p — l\fp— 1 + k 

k 



(47) 



^ fc (2p-l)(2p-2)--.(2p-fc)p 2 (p+l) 2 ---(p + fc-l) 

^ l J fc! A;! 2 

fe=i 



^ 2 V 7 / V 7 

fc=l j=i v J 7 i=l v J 



Since 

k 



f[ (l - ^ = 1 - 2j9# fc + 4p 2 #({l} 2 ; fc) = 1 - 2p# fc + 2p 2 (# 2 - #f ) (mod p 3 ) 
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and 

fc-l , N 2 fc-l 



n f 1 + -) = n f 1 + - + - 1 + + p 2H i 2 -i + v^({i} 2 ; * - 1) 

7=1 ^ 3 ' 7=1 ^ "7 / 



= 1 + 2pH k ^ + 2p 2 if 2 _ 1 - p 2 H { k 2 \ (mod p 3 ) , 

substituting these congruences in (j47p and simplifying we obtain 
(48) 

Note that 

p-1 rr(2) p-1 fc-l p-1 fc 

fc=i fc=i j=i J fc=i i=i J 

i=i J j=i J 



and therefore 

H k-i 



p- 1 tr( 2 ) . » 2 



(49) 2 ^^1 = ( H «) 



fc=l 

Now by Lemma [31 from (HH1) and (H9l) we have 

/,_J x 2p - 1\ /2p - k - ? x 2 

p + k J \ p — 1 — k 



(50) E(-D*(*. ,!V 2 f * f] ^ItW (mod/). 

fc=0 V ' 

From Lemma [5] we find easily 
(51) 



W =4(l + 2pi/ p _ 1 -j9 2 ^ 2) 1 ) (modp 5 ). 



Now from Pol) . (|50|) and (151~|) . by Lemma [3j for any prime p > 5, we have 

]T(-l) fc W 5 (205A; 2 + 160A; + 32) = 32p 2 (l + p 2 H^\) (l + 2pH p _ 1 - p 2 H^ 

fc=0 ^ ' 

= 32p 2 + 6Ap 3 H p ^ 1 (modp 7 ). 

The validity of this congruence for p = 3 can be easily checked by straightforward 
verification. Taking into account that for an odd prime p, 

2k\ 5 

J =0 (mod p 5 ) 

for k = ^rr, . . . ,p — 1, and applying Lemma [3j we get the second congruence of Theo- 
rem [HI □ 
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6. A SUPERCONGRUENCE ARISING FROM A SERIES FOR THE CONSTANT K. 

In Section [3] we proved the accelerated formula 

;-27) n - 1 (15n-4) 



n=l 



n 



3 ( 2 :) 2 ffl 



which was earlier conjectured by Z.-W. Sun. Motivated by this series, Z.-W. Sun 
Conj. 5.6] formulated the following conjecture on supercongruences: for any prime p > 3 
and a positive integer a, 

P ^ I5k + 4f2k\ 2 (3k\ Jp a \ a . J 2+a . 

In this connection we prove here the following theorem. 
Theorem 9. Let p be a prime greater than 3. Then 

Proof. Consider the WZ pair {F,G) defined in the proof of Theorem [1] with a = 1/3, 
f3 = 2/3, a = : 

_ (-l)"(3n + l)!n! 3 (l)l(l) 2 k (n + 2k + l) 
{,) (2n)!27" ' m 2 m 2 

V3/n+fc+l V3/n+fc+l 

„, ,* „, 45n 2 + 63n + 22 + 18A;(3n + fc + 2) 

Gin, A; = Fin, k) ■ ; — — . 

V ' ; V ' ; 18(2n + l)(n + 1 + 2k) 

Then it is readily seen that the pair (F, G) given by 

— . , N (-27)- fc (2A;)!(3£; + 3n)! 2 n! 2 , JW , 

« = (3t+l)ltP(t + «)P(3n)P « 15fc+4 )( 3fc + 1 ) + 1M3* +»+D), 

,, 3A; 3 (3A;- l)(-27) 1 - fc (2A;)!(3fc + 3n)! 2 n! 2 (2n + A; + 1) 
G(n, k) = 



(3£;)!£;! 3 (fc + n)! 2 (3n + 2)! 2 
is its dual WZ pair (see Proposition C or [20, Ch. 7]), for which we have 

(52) F(n + l,k) - F(n, k) = G(n, A; + 1) - G(n, k), n, k > 0. 

Summing ( 152]) over k — 0, 1, . . . , p — 1 and observing that G(n, 0) = 0, we obtain 

p— i p— i 

(53) ^ F(n + 1, fc) - J2 Hn, k) = G{n, p). 

k=Q k=0 

Further, for every integer n satisfying < n < we have 

(54) Gin p) = (-27) 1 -n2p)!(3 P + 3n)! 2 n! 2 (2n + P+ l) 

1 J 1 (3p-2)!(j9-l)! 2 j9!(p + n)!2(3n + 2)!2 " U lm ° a P J ' 
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since the numerator is divisible by p 8 and the denominator is divisible by p 5 and not 
divisible by p 6 . Now from (j55|) and for any non- negative integer n < we have 

p— 1 p— 1 p— i 

(55) ^F(O,A0 = ^F(l,A0 = "- = ^F(n + l,Jfe) (modp 3 ). 

fc=0 k=0 k=0 

Moreover, from ( |55|) we obtain 

Efcl^fiM) (modp 3 ), if p = l (mod 3); 
Efc^fir^) (modp 3 ), if p = 2 (mod 3). 



(56) F^fc) 
Since 



?*M>-?^(?)"(? 



fc=0 fc=0 

it is sufficient to prove the required congruence for the right-hand side of f|56|) . We 
consider separately two cases depending on the sign of (|). 
First, let p = 2 (mod 3), then we have 

<p-2 ^ _ <-»)-W(3* + p-2)Pfi*)P 



3 ' / (3fc + l)!fc! 3 (fc+2^)! 2 (p-2)! 2 



Note that if 1 < A; < ^r^, then the denominator of F(—^-, k) is not divisible by p and 
the numerator is divisible by p 2 . Therefore, we have 

P-2 ,A _ n _ 2 x , 1 n P-2 



(57) F^i-^— ,fcj=0 (modp 2 ), fc=l,2, 

Let ord p n be the p-adic order of n that is the exponent of the highest power of p 
dividing n. It is clear that 

(58) ord p(EzlP±l) = ord 3^mK2p-l)P(¥)'' _ , 

Similarly, considering the disjointed intervals < fc < < k < , and 

< A; < p — 1, we obtain 

(59) ord p F ( ^^,A; ] > 3, fc = > + X > • • • >P ~ L 



3 

Thus from f l57|) -f l59|) we have 
p-i 



fc=0 



or 



— /p-2 \_ 3 -^ 1 (^)!(2p-l)! 2 (£^)! 2 



/ — 2 \ 

= 4 + 2(p-2)(p+l) + 

fc=0 ^ ' 



(p + 2 )!(£±i)!3 (^i)!2 (p _ 2 )!2 

„2N 



x ((5p + 9)(p + 2) + 8(p - 2)(p + 1)) (mod p 2 ). 
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Taking into account that 

(2p-l)! 2 (p-l) 2 p ip+ l) 2 (p + 2) 2 ---(2p- If 



(p-2)! 2 (p + 2)! (p + i)(p + 2) (p-l) ] . 



1 '' (p — 1)! (mod p 2 ) = — (modp 2 ), 



(p+l)(p + 2) ^ ' ^ - ' 2 

(5p + 9)(p + 2) + 8(p - 2)(p + 1) = 2 (mod p), 
and simplifying, we get 

^ n /p-2 \ 2-3~ p -p! , , „, 

For primes p > 3, by Fermat's theorem, we have 3 P = 3 (mod p) and therefore, 
' ,_J X P-2 ,\ _ n 2 p! 



fc=0 

Now put n = ^i- and note that 



(60) E f (V' fc ) E ^ + 3 (y)! ( ^)i (m ° dp2) - 



Pj _ p(p-l)---(p-7l+l) (_l)n-l / p 

(61) (¥)!(f)! W n! P n ±jl 7 

~~ -p (mod p 2 ) = — 3p (modp 2 ), 



n 

where in the last equality we used the fact that n — 1 = is odd (otherwise, a prime 
p > 3 must be even). Finally, substituting ( 16T|) in ( 160|) we obtain 

p 1 / 2 \ 

Y^F\X^-,k\ =-2p-2p=-4p = 4(|jp (modp 2 ), 

fc=0 ^ ' 

which proves the theorem for primes p = 2 (mod 3). 
Now suppose that p = 1 (mod 3). Then we have 

]T F(0, k)^F (tzl, k) (mod p 3 ), 



fc=0 fc=0 

where 

p-1 \ _ (-27)- fc (2A;)!(3A; + p-l)! 2 (£3i)! 2 



(3A; + l)!A;! 3 (A;+£^l)! 2 (p-l)! 2 



((15fc+4)(3fc+l)+2(p-l)(9fc+p+2)). 



Note that if 1 < k < then the denominator of F(^-,k) is not divisible by p. On 
the other hand, ord p (3k + p — l)! 2 = 2 and therefore, F(^^,k) = (modp 2 ). It is 
clear that 



1 3 3 
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Similarly, if £±2 < k < &=1, then ord p F(^, k) > 3, and if < k < p - 1, then 
ord p F(^i, k) = 3. Therefore, we have 

(62) 'f/^^F^^F^/-^) (mod A 
where 

(63) f(^^,0] = 4 + 2(p- l)(p + 2) = 2p (modp 2 ) 



and 



F 



3 

p-1 p-l^ 3 1 - p (2p-2)! 2 (p(5p- 1) +2(p- l)(4p- 1)) 

(p — 2)!=p! (modp 2 ), 



Noting that 

(2p-2)! 2 p (p+ l) 2 (p + 2) 2 ---(2p-2) 2 = p 
p!(p— l)! 2 p— 1 (p — 2)! p — 

and applying Fermat's theorem, we have 

p — lp — 1\ _ 2 • p! 



3 ' 3 ;-(^)!(^)T {m ° dp2) - 
Now setting n = £jp and taking into account that 

we have 

(64) F ( ) = 2p (mod p 2 ). 

Now by (IS2"|l - (IS4"jl . we obtain 

p-i 



5=1 



£jr =4p = 4(|)p (modp 2 ) 

fc=0 ^ ' 



which proves the theorem for primes p = 1 (mod 3). □ 
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